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ry^ ' Context. Rapidly rotating stars show short-period oscillations in magnetic activity and polar appearance of starspots. 

. , Aims. The aim of this paper is to study large-scale shallow water waves in the tachoclines of rapidly rotating stars and their connection 

r^^ . to the periodicity and the formation of starspots at high latitudes. 

(Ph' Methods. Shallow-water magnetohydrodynamic equations were used to study the dynamics of large-scale waves at the rapidly rotating 

stellar tachoclines in the presence of toroidal magnetic field. Dispersion relations and latitudinal distribution of wave modes were 
derived. 

Results. We found that low-frequency magnetic Rossby waves tend to be located at poles, but high-frequency magnetic Poincare 
~^ ' waves are concentrated near the equator in rapidly rotating stars. These results have important implications for the evolution of the 

stellar wind in young Sun-like stars. 

Conclusions. Unstable magnetic Rossby waves may lead to the local enhancement of magnetic flux at high latitudes of tachoclines in 
rapidly rotating stars. The enhanced magnetic flux may rise upwards owing to the magnetic buoyancy in the form of tubes and appear 
^ ' as starspots at polar regions. Magnetic Rossby waves may also cause observed short-term periodicity in the stellar magnetic activity. 

JL. ' Key words. Stars: magnetic field - Stars: activity - Stars: starspots - Stars: oscillations - Physical data and processes: 

^1^ ' Magnetohydrodynamics (MHD) 

rn 

C^ 1. Introduction 

o : 

Rotation and convection lead to stellar dynamo activity, which results in the strong concentration of magnetic fields on the surface. 

These magnetic features are called starspots, and they represent the sunspot analog for other stars. Interested readers may find 

detailed information about starspots in two recent comprehensive reviews: Berdyugina (I2005I I and Strassmeier (12009b . Hall (119721) 

• »— I explicitly postulated the starspot model as explaining the wave-like features in the light curves of rotating RS CVn stars. Doppler 

^^ imaging technique showed spots on the RS CVn star HR 1099, and they exhibit a similarity in shape and location to X-ray images 

H ! of solar coronal holes (Vogt and Penrod. ri983l l. Later, Strassmeier (1 19901 1 observed a cool polar spot on the rapidly rotating G5 

. . . ■ subgiant in the single-lined RS CVn binary HD 26337. The spot was rather long-lived, persisting for at least the nine years of 

observations. 

Zeeman-Doppler imaging of the circularly polarized spectrum has led to the first magnetic maps of the active Kl subgiant 
in the RS CVn binary system HR 1099 = V711 Tau (K1IV-I-G5V) other than the Sun, which also showed the polar location of 
starspots (Donati et al. 119921 ). It was generally accepted that starspots on magnetically active, cool stars mainly appear near the 
poles. Schiissler and Solanki (|1992| l suggested that the preference for high latitudes is due to the rapid rotation to these stars, which 
in turn leads to a dominance of the Coriolis force over the buoyancy force in the dynamics of magnetic flux tubes. Later, Schiissler 
et al. dl 9961 1 caiTied out simulations on solar-type stars in an effort to predict the emergence latitude of magnetic flux tubes, after 
evolution from a toroidal magnetic field located at the base of the convection zone. They carried out simulations for stars of different 
ages and equatorial rotation velocities. The results predict high-latitude spots since, in a fast rotator, Coriolis forces acting on a 
rising magnetic flux tube cause movement parallel to the rotation axis. 

On the other hand, Barnes et al. (il998) studied two rapidly rotating G dwarfs in the a Persei cluster, from spectra taken at 
the William Herschel Telescope. The image reconstructions demonstrate that both these stars have cool high-latitude regions or 
polar crowns and low-latitude features, contradicting the suggestion that only high-latitude spots should be present. The magnetic 
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topology of AB Dor (young KO dwarf) is also found to be very complex, with at least 12 different radial field regions of opposite 
polarities located all around the star (Donati et al. 1 19991 [20031 1. Significant azimuthal field fluxes are also detected in the form of one 
negative polarity region close to the equator, a series of positive polarity patches at intermediate latitudes and an almost complete 
ring of negative polarity encircling the rotational pole at high latitudes. Starspots on the young solar analog EK Dra (rotation period 
of 2.68 days) are observed to be formed everywhere, from the low latitudes up to polar regions (Strassmeier 2009). Therefore, the 
current view is that starspots appear at all latitudes in rapidly rotating stars. 

The high-latitude sunspots could change the structure of coronal magnetic field as well as the formation of prominences, which 
may lead to the predominance of solar flares and coronal mass ejections (CMEs) near the polar regions. CMEs are believed to have 
tremendous influence on the atmosphere of exoplanets, leading to their erosion (Khodachenko et al. 120071 ). The polar CMEs may 
propagate significantly out off the ecliptic plane, which may cause more convenient conditions of planetary evolution; therefore, the 
detailed study of starspot distribution in rapidly rotating young solar analogs is very important for planetary evolution as well. 

Observation of starspots at all latitudes complicates any explanation of their appearance. There are at least two suggested 
mechanisms of starspot formation in rapidly rotating stars. The first mechanism, which we already mentioned above, is that the 
magnetic tubes tend to move to the polar regions owing the Coriolis force acting during the rising across the convection zone 
(Schussler and Solanki [19921 Schussler et al. [19961 DeLuca et al. [19971 Granzer et al. [20001 Holzwarth et al. [20061 Isik et al. [2007l l. 
The second mechanism is the transport of photospheric magnetic flux towards the polar regions by meridional advection (Schrijver 
and Title I200TI Mackay et al. 12004b . The first mechanism may explain the starspot appearance at high latitudes, but it may fail to 
explain the equatorial spots because the magnetic tubes should all tend to higher latitudes in rapid rotators. On the other hand, the 
second mechanism may explain the distribution of spots at all latitudes, but it requires a strong meridional flow in order to carry the 
magnetic flux towards the poles. 

Another interesting observational indication is the cyclic behavior of stellar activity with shorter periods than the solar cycle (~ 
1 1 yr). Massi et al. ([1998 2005) report short-term activity cycles on UX Arietis, which is a binary system formed by a spotted KO 
IV primary and a G5 V secondary with an orbital period of 6.44 days. Short-period activity cycles have been observed on r Bootis 
(Donati et al. I2008I I. Olah et al. (I2009I I analyzed 20 active stars and came to the conclusion that stellar activity cycles are generally 
multiple and variable. Recently, Lanza et al. (I2009I I have observed a short-term cyclic oscillation of the total spotted area on the 
photosphere of CoRoT-Exo-2a, which is supposedly caused by Rossby waves. The rotation period of CoRoT-Exo-2a is 4.522 + 
0.024 days, while the total spotted area shows a cyclic oscillation with a period of 28.9 ±4.3 days. Almost all observed periodicities 
are shorter than the solar cycles and are several times longer than the stellar rotation periods. It should be mentioned, however, that 
solar activity also undergoes intermediate periodicities in the range of several months (Rieger et al. 11984! Carbonell et al. 11990! 
Oliver et al. '1998) and ~ 2 years (Sakurai [1981! Gigolashvili et al. 11995! Vecchio et al. 1201 Oi l. These periodicities are explained 
in terms of magnetic Rossby wave instability in the solar tachocline that is the joint effect of latitudinal differential rotation and 
toroidal magnetic field (Zaqarashvili et al. l2010all20I0bl l. The coincidence between the solar short term oscillations and the cyclic 
behavior of spotted stars indicate that magnetic Rossby waves can be important in rapidly rotating stars. 

Here we study the behavior of magnetic shallow water waves in tachoclines of rapidly rotating stars and their role in the 
latitudinal distribution of starspots. The tachocline is the thin layer located between the external convection envelope and the internal 
radiative core (in the case of the Sun see Spiegel and Zahn !1992l l. Rossby waves are well studied in the geophysical context (Gill 
11982b : however, the presence of magnetic fields significantly modifies their dynamics (Zaqarashvili et al. 12007! !2009! Heng and 
Spitkovsky !2009! Mathis and de Brve l20lTb . We have used the shallow water magnetohydrodynamic (MHD) equations (Oilman 
12000b to study the latitudinal structure of magnetic shallow water waves in tachoclines of rapidly rotating stars. We first write the 
governing equations of shallow water MHD in the rotating frame. Then, we derive the dispersion relations and latitudinal structure 
of magnetic shallow water waves in different configurations of toroidal magnetic field, which is followed by the discussion of the 
connection between the shallow water waves and the starspots at stellar photospheres. 

Finally we discuss the implications of our findings in the context of the evolution of young stellar winds of Sun-like stars. There 
is observational evidence that zero-age-main-sequence solar-type stars rotate more than ten times faster than the present Sun (e.g., 
Giidel et al. 11997! Ribas et al. !2005! Lammer et al. 1201 lb . As a consequence, young solar-type GV stars, including the early Sun, 
should have vigorous magnetic dynamos and correspondingly strong high-energy emissions. From the study of solar-type stars with 
different ages, Skumanich (!1972b , Simon et al. (!1985b and others have shown that these young stars lose angular momentum with 
time via magnetized stellar winds, thus leading to a secular increase in their rotation period (Durnev ri972b . This rotation slowdown 
is well fitted by a power law roughly proportional to f'^^, where f is the stellar age in units of Gyr(e.g., Skumanich[T972 ; Soderblom 
!1982! Avres [T997l i. As a response to slower rotation the solar dynamo strength diminishes with time, causing the Sun's high-energy 
emissions also to undergo a significant decrease (Ribas et al. 12005b . The fundamental question of whether the young fast-rotating 
Sun and similar stars have experienced some periods of stronger magnetic activity has a major impact on the evolution and plasma 
interaction of planetary atmospheres in the Solar System and beyond. 

2. Governing equations 

We use a spherical coordinate system (r, 9, <p) rotating with the stellar equator, where r is the radial coordinate, 6 is the co-latitude, 
and <p the longitude. We consider an unperturbed toroidal magnetic field B^ - B(8) sin 8. Then, the linearized form of shallow-water 
MHD equations can be rewritten in the rotating frame as (Zaqarashvili et al. 12007! [20 lOab 

2i2() cos 011,1, H sin9 h 2 cos 9bd, = 0, ( 1 ) 

dt ^ Ro 89 AnpRa d(f> AnpRo * 

du^ oo ^-^ S dh bBsm9dB{9) B(9) db^ 2B(9) . ^ 

— — -I- 2i2n COS 9ug + COS 9bg = 0, (2) 

dt Ro d(f> A-npRo 89 AnpRo d(f> 4-KpRo ^ ' 
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. 2 „dh Ho . due Hq du^ 

sin 6 — + — sin + = 0, 

dt Rq 86 Ro dcf> 

dbg B(e) diig ^ 



dt Ro d(p 



3 

(3) 

(4) 



db^ 1 d 

dt Rq dO 



(5) 



where ug - sin 0ug, ii^ - sin 6u^, bg - sin Obg, and b^ - sin 6b^ are the velocity and magnetic field perturbations, Hq is the tachocline 
thickness, h the perturbation of the thickness, p the density, Rq the distance from the stellar center to the tachocline, g the reduced 
gravity (due to subadiabatic temperature gradient; Gill 1 19821 ) and Q.q the equatorial angular velocity. For a zero magnetic field, this 
system transforms into the HD shallow-water equations. 

We now perform a Fourier analysis of the form exp(-/(yf + is(p) and define 



= ^, = e, vJu.), mu), COS0 - u, -sin0-— - (1 - u )-— - D. 

2Qo gHo Anp ''^^'' 2Q.qRq ^'' ^ dO ^ ^ ' dfi 

After some algebra, we get a single equation (Zaqarashvili et al. 120091 ): 



(6) 



{AD + sn) 



s--eAHl-ii2) 



AD- s^- —s^A{D + 2^i) 
A~ 



iig - {A~ - jj. )ug + s a iig + 2a /uDug + jjs — (D + 2^)ug + iiD{a )ug - 0(7) 

A 



This is the general equation for wave dynamics in a shallow-water MHD system. It contains various kind of waves (magneto-gravity, 
Rossby, Kelvin). 

Equation |2]can be rewritten as 



eA^ + 2ea^sA 



S + 2A 



-a Ae/i + 



2A^en 



2as da 



M^ 



^2 _ Q,2^2 ^ ^2 _ ^^2 (1 _ ^2) ^2 _ ^,2^2 ^^ 

2eA^s+4ea^s^A^ 



I IN dug 



2;i 



T^ 



A^-a^s^ '^ {A2-a^^s'^)\_s^-eA2{\-^i2)Y A^-a^s^d/j 
Now substitution of a new function 



-, 2aeA^ da i 9\ 25^0- da 



Ug - Ugi exp 



\s 



S + 2A 



A^ — a^s^ 



a Aefi + — 



2A^en 



2as da 



52-6^2(1-1^2) A^-a^s^dfi 



dfi 



A- - a^s^ d^i 



1-^2 
tig = 0, 



9 sA + 2a s 
+ eA' - 



X"- - Q'2i2 



(8) 



(9) 



gives 



d (^ 2\ dugx 



l-yL/2 



+ eA'- 



sA + 2a^s^ eA^ + 2ea^sA 



2eA^s + Aea^s^A^ 



2 S + 2A 2 1 2 



A^ — a^s^ 



A^-a^S^ ^ "^ (^2_Q,2^2)[^2_e^2(l_^2)]A' ^2_„2^2 



T^ - 



ea A/i - 



2eA^ 



y + 



exl2 



s2-e^2(l_^2)'- i2_^^2(l_^2) 



/ -,\ I s + 2A 2 



2 A~ ~ a^s 



2 ,.2 



-'('-''')-i7l^^«'-<V(i-M=)- 



(s + 2A)e^a^A 



2^,2 ,3 



(^2-a2^2)[^2_f^2(l_^2)] 



^"(^"^')~ 



3e2^4 



[s^ - eAHl - n^)] 



^//-(l -^r) 



4 (^2 _ Q.2j2)2 

2£aA^s^ 



eA-4s da i j\ s'A}- 

+ -. — —sa—^l{l-^l} 



a'c-''^) 



as^ d^a I 2\ 



(^2 - a,2^2)[^2 _ g^2(l -^2)]^ 5^ 
Ug\ = 0. 



5^^^(l-^') + 



(10) 



i2-ff2^2 Q^t-\ ^> (^2 _ 0,2^2)2 \5ju/ ^ "^f A^-a^s^dfi^ 

At this stage we should fix the magnetic field profile. We consider two different cases of the magnetic field B^ - Bo -\Jl - jj~ 

and B^ = Bq/j ^fl^^Jp■, where Bq is constant parameter Here, B^ - Bq -\J1 - jj.^ means that the resulting magnetic field is at its 
maximum at the equator and tends to zero at the poles. The magnetic field is opposite in the northern and southern hemispheres. 
Then, B^ = Bq/i -\J1 - ii~ means that the magnetic field tends to zero at the solar equator and poles, reaching the maximum value 
at middle latitudes. This configuration is similar to the structure of the solar magnetic field, where the sunspots appear at middle 
latitudes (Oilman and Fox 1 19971 ). We now consider both cases separately. 
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Fig. 1. Left: normalized frequency A = a)/2Q.o vs ratio of Alfven to surface speed a = va/2Q.oR(). Here s = 1 (toroidal wave number) 
and £ = AQ^R^I gH{) - 100 (ratio of surface to internal gravity speed). The blue, red, and green asterisks (circles) correspond to 
the positive (negative) solutions of the dispersion relation Eq. |23]for the v = 0, v = 1, and v - 2 harmonics, respectively. Right: 
normalized frequency A - (jj/2Qo vs e for 5=1 and a-0. 1 . 



3. Uniform magnetic field {B^ - Bq ^Jl - fi^) 

The equation B^ = Bo ■\Ji - ^J^ implies a - const and Eq.[TO]is rewritten as 

+ Bi+ B2i? +£3(1- fJ^) + B41? (1 - 1?) Ugi = 0, 



|;('-''')t^ 



l-jj2 



where 



sA + 2a'^s^ 



Bi = eA^ - 

eA^ + 2ea^sA 

Bi -^ 7^— + 

A- - a-s-^ 

eA^ 



2eA^s + Aea^s^A^ 



2A 



(/12 - a^-s^)[s^ - eA\\ - ^i^)] A^ - a^s^ 
1 S + 2A 



ea'-A- 



2eA^ 



s^-eAHl-i^^y 



H 



-ea'^A, 



s2- 6^2(1 -^2) 2^2_Q,2j2 

2 (s + 2A)e^a~A^ 



1 is + 2A) 2 4,2 



le'-A 



2i4 



(^2 - ff2 j2)[^2 _ ^^2(1 _ ^2)] [j2 _ ^^2(1 _ ^2)]2 " 



(11) 

(12) 
(13) 
(14) 
(15) 



We now consider the approximation of rapidly rotating stars or a spherical shell with reduced gravity (such as in the overshoot 
region): 



eJ^^»X. 



gHo 



(16) 



In this case Eq.[TT]is a spheroidal wave equation, but with a slightly compUcated form (for example, the terms with yu"* arise). An 
analytical method for the solution of such equations was proposed by Longuet-Higgins (II 96511 968l l for the Earth's atmosphere. We 
closely follow this method. The equations contain two different branches of solutions: high and low frequency ones. The higher 
frequency branch is represented by magnetic Poincare waves and the lower frequency one by magnetic Rossby waves. We study 
both cases separately. 



3.1. High-frequency waves 

When e » 1 , Eq. [TT]is most easily satisfied if A-, a^, and fr are small enough to compensate for the large e. Let us introduce a new 
independent variable ;; = e'^^ju and consider that A ~ e '''^L and a ~ e '^^L, where L is near unity. Keeping only the terms > e'^"^ 
in Eq.[TTl we get 

(17) 



(18) 



\d^ A ,1 


Mfll -- 


where 


--'44 


ea s 
A 
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Fig. 2. Normalized frequency A = lij/2D.o for four different values of e (10, 100, 1000, 10000). Here e is shown on a logarithmic 
scale. The blue, red, and green asterisks correspond to the positive solutions of the dispersion relation Eq. |23]for the v = 0, v = 1, 
and V = 2 harmonics (with s = 1), respectively. The dashed line corresponds to A = e"'^^ slope. Here the ratio of Alfven to surface 
speed a = va/2Q()R^) equals to e"'^- as assumed in the derivation of Eq. [TTl 



s=1,a=0.1,E=100, A=0.72012 



1,a=0.1,E=100, X=0.57518 




s=l,a=0.1,e=100, ?i=0.36686 




Fig. 3. Numerical solutions of Eq. [TT]for three different values of the normalized frequency A = w/2Qo (solid lines), jj - cosO, 
where is the co-latitude, hence /y = corresponds to the equator and j-i = I corresponds to the northern pole. Upper, middle and 
lower panels correspond to the harmonics with v = 2 (/1=0.72012), v = 1 (/1=0.57518) and v = (/l=0. 36686), respectively. Dashed 
lines represent corresponding analytical solutions from Eq. |20] All solutions are concentrated at low latitudes. We take .? = 1, 
£ = 100 and a = 0.1 during the simulations. 

This is the Weber (parabolic cylinder) equation, and the finite solutions atrj ^ +co yield (Abramowitz and Stegun [r964l) 



-- =2v+l (v = 0,1,2,...) 
and then 

llgi cc e^'-'^' Hy{T]), 

where //,,(/;) is the Hermite polynomial v. The function uo[ is exponentially small beyond the turning points 



T]^ V2v+ 1, 
that is to say. 



= e-'/4V27TT. 



(19) 



(20) 



(21) 



(22) 



Due to the large e, the co-latitude satisfying this condition is close to 90*^, therefore the solution is confined to the neighborhood 
of the equator 

Equation [T9]yields the wave dispersion relation as 



^-:" 



£a 



0. 



(23) 



For non-magnetic case, this equation transforms into Eq. (8.11) of Longuet-Higgins (I1968I I. 
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For high values of e, the solutions of this cubic equation can be given as 

(2v+ 1)1/2 ^ I ^^2X 

l-^l (25) 



(2v+l)ei/2\ 2 

The first solution represents the high-frequency branch of shallow-water equations, i.e. magnetic Poincare waves (or magneto- 
gravity waves). For a -Q, expression (l24l i is the same as for the nonmagnetic case (Longuet-Higgins |l968l ). The dispersion relation 
yields A ~ e '/"^, so it agrees with the approximation used in the derivation. On the other hand, expression ( |25] ) represents the 
magnetic Rossby waves, but it does not agree with the approximations as it yields A ~ e '^2. Therefore, the second solution does 
not describe the wave dynamics correctly. 

We solve the dispersion relation Eq.|23] numerically for different values of the ratio of Alfven to surface speed a = va/2Q.oRo 
and the ratio of surface to internal gravity speed e = 4-Q.^R^/gH(). Stability analysis of shallow-water MHD equations in the presence 
of latitudinal differential rotation shows that the s = I harmonic is the most unstable one in all considered cases (Oilman and Fox 
119971 Zaqarashvili et al. l2010al2010bl) . Therefore, in what follows, we always consider the 5=1 harmonics. 

The normalized frequency A = a)/2£lo vs a and e is shown on Fig. 1. The left hand panel shows the dependence of /I on a for 
s - 1 and e - 100. We see that the frequency depends only slightly on a, i.e. on magnetic field strength. The right hand panel shows 
the dependence of /I on e for s = I and a = 0.1. The dependence of the wave frequency on e is more pronounced; higher e yields 
smaller A. The figure displays only those two solutions of the dispersion relation Eq. |23]that satisfy the condition A ~ e ''"*. The 
third solution, which has lower frequency and does not satisfy A ~ e '^"*, is not shown here. 

The normalized frequency A for four different values of e (10, 100, 1000, 10000) is shown on Fig. 2. The figure shows that the 
increase of e leads to a significant decrease in the normalized frequency. For example, e = 10 yields Tosc ~ 0.8 Trot, while e - 10000 
yields Tg^^ ~ 5 T^ot- Here Tqsc is the oscillation period of s = 1, v = harmonics, and Tiot is the stellar rotational period. This means 
that faster rotation leads to a longer oscillation period than the rotational period. 

To check the analytical solutions, we solved Eq.[TT]numerically for several frequencies, which were found from the dispersion 
relation Eq.[23] We searched the bounded solutions of Eq.[TT] varying A close to its analytical counterpart. The numerical solutions 
of Eq.[TT]for three different values of A are shown on Fig. 3. We see that all solutions are bounded and concentrated at low latitudes 
as suggested by the analytical solution. The frequencies for which the bounded solutions are calculated are in good coincidence with 
those obtained from the analytical dispersion relation Eq.|23] Therefore, the wave modes with frequency A ~ e^^^* are concentrated 
at low latitudes. These modes are magnetic Poincare waves (Zaqarashvili et al. l2007l l. 

3.2. Low-frequency waves 

We now concentrate on lower frequency waves by considering A - e '^^L and a - e '^^L, where L is close to unity. Then eA^ ~ 1, 
which means that s^ - 6/1^(1 - fi^) tends to zero for s - I and small fi. As a result, some terms in Eq.[TT]tend to infinity. Therefore, 
Eq.[TT|is most easily satisfied if 1 - //^ is small instead of /i^. This means that the solutions are located near the poles not near the 
equator. We again follow the steps of Longuet-Higgins. Returning back to the variable 6, Eq.[TT]is written 

sin 0— (sin 0-^1 + [-s^ + Bi sin" + B2 cos^ flsin^ + B3 sin'^ + B4 sin"* Scos^ 0]ugi = 0. (26) 

00 \ 00 



For small we have 

d^Ugi 1 dugi 

56*2 '^ 80 '^ 



^ + Bi + B2(l - 0~) + Bi0^ + Ba0\1 - 0^) 



«fli = 0. (27) 



In order to estimate different terms in this equation, let us consider - e '^^77, where t] is near unity. Keeping only the terms 
> e''2 we get from Eq.|27] 

(28) 



(29) 
(30) 

(31) 



d-ue\ 1 dugx 
80- ^ 80 ^ 


-^+A,-A202 


Ugi = 0, 


where 


\ '^ 


6/^2 + 3ea~sA 

^2 - ^2^2 + 
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Now if we consider mode frequencies in the range 
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Fig. 4. Left: normalized frequency A = ai/lQo vs ratio of Alfven to surface speed a = va/2Q.oR(). Here s = 1 (toroidal wave number) 
and e = AQ^R^I gH() = 100 (ratio of surface to internal gravity speed). The blue, red, green, and yellow asterisks (circles) correspond 
to the positive (negative) solutions of the dispersion relation Eq. [36]for the n - 1, n - 2, n - 3, and n - 4 harmonics, respectively. 
Right; normalized frequency A - w/2Qo vs e for s = I and a = 0.1. 
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relation Eq.[36]for the « = l,n - 2, n = 3, and « = 4 harmonics, respectively. Dashed lines correspond to A = e '^^ slope. Here the 
ratio of Alfven to surface speed a = va/2Q.^)Ro equals e"'^^ as assumed in the derivation of Eq.l27l 



but A + ±as, then Ai is positive, and we can define a new variable x - V^a^^, which gives the equation 



x——r- + -r: — + 



X 2n + s + 1 s 
~4 "^ 2 4l 



dx^ dx 
where n is a nonnegative integer, and 

2{2n + s+\) 



ue\ - 0, 



(32) 



(33) 



This equation is self-adjoint and its solutions are sometimes called Laguerre functions ijj'ljix) (see for example, Arfken and Weber 
12005b . These functions are bounded at x = and x — oo and are linked to the associated Laguerre polynomials as 



uei = uei{Q)x'-''e '-'mx), 
where 

Zj {s + k)\{n - k)\k\ 

k—{) 

so that i > - 1 . 

Then Eq.[33]gives the wave dispersion relation as 
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Fig. 6. Upper panels: Numerical solutions of Eq. [TT]for two different values of e and corresponding normalized frequency A = 
w/2Qo. ju - cos 6, where 6 is co-latitude, hence // = corresponds to the equator and j-i = I corresponds to the northern pole. We 
take s - 1 and a - yfe during the simulations. Lower panels: Corresponding analytical solutions plotted from Eq. [34l left panel 
shows the solution with n - \ and the right panel shows the solution with « = 3. 



The sign of A - as and A define the sign of the righthand side of Eq.| 



A^ — a^s^ 



o(:i^<o)' 



yields a + (-) sign. We solve the dispersion relation Eq.[36]numerically for different values of a and e. 

The normalized frequency A\s a and e is shown on Fig. 4. The left hand panel shows the dependence of /I on a for s - \ and 
e = 100. We see that unlike the previous case here the frequency significantly depends on a, i.e., on magnetic field strength: higher 
a yields larger A. The right hand panel shows the dependence of /l on e for s - \ and a - 0.1. The dependence of wave frequency 
on e is less pronounced here. 

The normalized frequency A for four different values of e (10, 100, 1000, 10000) is displaced on Fig. 5 . The figure shows that 
the increase of e leads to a significant decrease in the normalized frequency. For example, e = 10 yields Tosc ~ l-75rrot, while 
e - 10000 yields Tosc ~ lOOOrrot. Here Tqsc is the oscillation period of the s - \,n - 1 harmonics and Trot is the stellar rotational 
period. This means that faster rotation leads to a longer oscillation period with respect to the rotational period. However, the increase 
of e leads to the deviation of A from e"'^^, while the analytical solution is obtained in the approximation when A is the same order 
of magnitude as e"'''^. 

To check the analytical solutions, we solved Eq.[TT]numerically for several frequencies, which were found from the dispersion 
relation Eq.|36] We searched the solutions of Eq.fTTIfixing A close to its analytical counterpart. The numerical solutions of Eq.fTTIfor 
two different values of e and corresponding normalized frequency A is displaced on Fig. 6 (upper panels). We see that all solutions are 
concentrated near the polar regions, as suggested by the analytical solution. The frequencies for which the solutions are calculated 
are in good coincidence with those obtained from the analytical dispersion relation Eq. [36] The lower panels of this figure display 
the corresponding analytical solutions plotted from Eq. |34] The figure shows that the analytical and numerical solutions show the 
same latitudinal behavior 

Therefore, the wave modes with low frequency, i.e. A ~ e '^^, are concentrated near the poles. These modes are magnetic Rossby 
waves (Zaqarashvili et al. l2007l l. 
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4. Non-uniform magnetic field (B^ - Bqh yj\- h^) 

In this case, Eq. [TO]is rewritten as 
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(37) 
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We follow the steps of the previous section and introduce a new independent variable 77 = e'^V' ^"d also consider A-e '^"^L 
and a - e^'^I^L, where L is near unity. Keeping only the terms > e'^"* in Eq.[32]we get 



where 



«fli = 0, 



2 2 
, s afX 



^2 



(45) 



(46) 



This equation and the corresponding dispersion relation was already studied in detail by Zaqarashvili et al. ("20091), so we do 
not repeat the calculations here. The solutions are concentrated near the equator as in the previous section. The resulting dispersion 
relation contains magnetic Poincare and magnetic Rossby waves. The influence of the magnetic field on Poincare waves is even less 
pronounced in this profile than in the previous section (compare Eq.[24]of the previous section to Eq. 11 in Zaqarashvili et al. l2009| l. 
The reason is that the magnetic field with this profile tends to zero near the equatorial regions, where the waves are concentrated. 
Therefore, it has almost no influence on the wave dynamics. 

However, we now concentrate on lower frequency waves considering A - e^^^^L and cq - e^^^^L, where L is near unity. Then, 
some terms in Eq.fTTItend to infinity for s - \ and small ^. Therefore, Eq.[TT]is most easily satisfied if 1 - ju2 is small instead of /i2. 
This means that the solutions are located near the poles in this case. Returning back to the variable 9 Eq.[32]is written 



sin 6*— I sin 6-^ \ + [-s^ + Bi sin^ 6* + B2 cos^ sin^ + B3 sin"* 6* + B4 sin"* Ocos^B + B^ sin^ 6» cos** + Be sin"* cos"^ 9+ 
09 \ 09 I 

+ Bj sin"* ecos'' 9]ugi = 0. (47) 

For small 9 we have 
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Mfli = 0. (48) 
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Fig. 7. Schematic picture of connection between Rossby waves and polar spots in rapidly rotating stars. Magnetic Rossby waves are 
concentrated near the polar regions in tachoclines of rapidly rotating stars. The erupted magnetic flux tubes lead to the appearance 
of polar spots at the stellar surfaces. 



To estimate different terms in this equation, let us consider 
we get from Eq.l48] 



e ^^'^T], where rj is close to unity. Keeping only the terms > e 
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(49) 



(50) 



(51) 



One can find that Eq. |49]is the same as Eq.l28] therefore the dispersion relations and solutions are the same as in the previous 
section. This means that the magnetic field structure has no significant influence on the dynamics of magnetic Rossby waves. This 
is probably because the difference between the two profiles of the magnetic field is most important at the equator but much less 
important near the poles, where the eigenfunctions of magnetic Rossby waves are concentrated. 



5. Discussion 

Magnetic field and activity determine the structure of stellar chromosphere/corona and have a significant influence on the evolution 
of exoplanetary atmospheres (Khodachenko et al. l2007l l. Therefore, to study the stellar activity/magnetic field is crucial for recently 
discovered exoplanets. Magnetic activity is connected to starspots, which appear as strong concentrations of magnetic flux at stellar 
surfaces and thus are similar to sunspots. Emerged magnetic flux is probably connected to deeper regions, namely to the tachocline, 
which is a thin transition layer between the external convection envelope and the internal radiative core. The tachocline may prevent 
the spread of the solar angular momentum from the convection zone to the interior (Spiegel and Zahn ll99"2l Gough and Mclntyre 
11998b . and it is probably the place where the large-scale magnetic field that governs the solar/stellar activity is generated/amplified. 
The magnetic flux rises upwards due to the magnetic buoyancy in the form of tubes and appears as spots at the surface. 

Starspots on rapidly rotating stars are observed to be formed everywhere, from low latitudes up to polar regions (Strassmeier 
I2009I I. Enhanced Coriolis force in rapid rotators may push the rising magnetic flux tubes, which begin to rise in lower latitudes as in 
the solar case, towards higher latitudes to form polar spots (Schiissler and Solanki [1992t Schiissler et al. |1996t DeLuca et al. 1 19971 
Granzer et al. I2000t Holzwarth et al. 120061 Isik et al. 12007b . But rapidly rotating stars show low latitude spots as well, which is 
difficult to explain as all tubes should tend to higher latitudes. A meridional flow may also transport the photospheric magnetic flux 
towards the poles (Schrijver and Title 2001), but the flow should be significantly enhanced in rapid rotators compared to the Sun 
(Mackayetal.|20Wll. 

Besides the longer periodicity caused by the a - ojor the a^ dynamo (~ 11 yr for the Sun), stellar magnetic activity also shows 
shorter periodicity (Massi et al. 119981 120051 Donati et al. 120081 Lanza et al. 120091 ). Almost all observed periodicities are much 
shorter than the solar cycles but longer than the stellar rotation periods. This fact fits well with Rossby-type wave activity. There are 
observations of intermediate periodicity in solar activity (Riegeret al. 119841 Carbonellet al. 1990i Gigolashvili et al. 119951 Oliver et 
al. 1 19981 Vecchio et al. 120 10b . which were explained by the instability of magnetic Rossby waves (Zaqarashvili et al. l2010all2010bb . 
The rapid rotation may enhance the role of large-scale Rossby waves in the dynamics of stellar interior and atmospheres. 
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Here we study the linear dynamics of shallow water waves in rapidly rotating stellar tachoclines with toroidal magnetic field. 
The latitudinal differential rotation is not taken into account because it significantly complicates the solutions (Dzhalilov & Staude 
120041 Mathis 2009, Zaqarashvili et al. 2010b). We derive analytical dispersion relations and find the dependence of the wave 
frequencies on various parameters, which are expressed by angular velocity, reduced gravity, and magnetic field strength. Then we 
plot the latitudinal structure of the first few harmonics. We find that high-frequency magnetic Poincare waves are located mainly at 
low latitudes (Section 3.1; Fig. 3). In contrast, low-frequency magnetic Rossby waves are located near the poles in rapid rotators 
(Section 3.2; Fig. 6). Heng and Spitkovsky ( I2009I I also find that the low-frequency waves, called magnetostrophic modes, reside near 
the poles. Actually, the magnetostrophic modes defined by Heng and Spitkovsky ( 120091 ) are the same wave modes as the magnetic 
Rossby waves defined here. The difference between these two considerations is the magnetic field structure. Heng and Spitkovsky 
( I2009I I used the radial magnetic field, which is valid for neutron star-related applications, while we use the toroidal magnetic field, 
which is applicable to stellar tachoclines. The latitudinal differential rotation and toroidal magnetic field may lead to the instability 
of magnetic Rossby waves, which may enhance the magnetic flux in the polar regions. Then the enhanced magnetic flux may rise 
upwards owing to the magnetic buoyancy and appear as polar spots (see Fig. 7). 

The important parameter governing the dynamics of shallow water waves is e and it should be estimated for rapidly rotating stars. 
The radius, density, and width of the stellar tachocline can be used from the solar values as Rq = 5 10'° cm, p - 0.2 g cm"^, and 
Hq = 0.01 Rq = 5 10*^ cm, respectively. The subadiabatic temperature gradient of the tachocline may significantly reduce the actual 
gravity. For example, the reduced gravity at the main, stably stratified part of the solar tachocline can be as low as 5 lO^cms"^ 
(Zaqarashvili et al. |2009t . The typical value of the equatorial angular velocity of the solar tachocline is ^20 = 2.7 lO^^s '. We 
assume that k is the ratio between angular velocities of a sample star and the Sun i.e. k = Qq/Qq- Then we have 

e= — ^-^x0.3k^. (52) 

Here k is significantly larger than one for rapidly rotating stars. For example, the young solar analog EK Dra (rotation period of 2.68 
days) yields /c a; 10, which leads to e a; 30. More rapidly rotating stars will lead to larger e. Therefore, the approximation of large 
e is applicable for most rapid rotators. It should be mentioned that the reduced gravity g and the tachocline thickness Hq are very 
sensitive parameters (there are no exact values even in the case of Sun) and thus may vary significantly from star to star. 

5.1. Polar spots, young stellar winds, and impact on early planetary environments 

As pointed out in the introduction, young solar-like stars rotate much more faster than the present Sun. With enhanced high-energy 
emissions and frequent flares, young stars are also expected to have more powerful stellar winds (e.g. Wood et al. 120021 Giidel 2007). 
Wood et al. ( I2001ll2002b have applied a method to infer stellar wind characteristics from observations between the interaction of the 
fully ionized coronal stellar winds and the partially ionized local interstellar medium. Modeling the associated Lyman-o- absorption 
features that are formed in the astrospheres of observed stellar sample stars have provided the first empirical estimate of coronal 
mass loss rates for GK main-sequence stars and indicated that the mass loss rates, so that the expected stellar winds appear to 
increase with stellar activity. 

However, more recent observations by Wood et al. (120051 ) of the Lyman-ff absorption signature of the astrosphere of about 0.55 
Gyr old solar-like star ^ Boo indicate that there may possibly be a high-activity cut-off regarding the mass loss in a radiation activity 
relation derived in the form of a power law by Wood et al. (2002). The observations of ^ Boo indicate that the mass loss of that 
particular star is about 20 times less than the average stellar mass loss value inferred from stars with an age of about 3.9 Gyr 

From the findings of our study we would expect a rotation-related modification of the outflow region for young stellar winds, 
which may have affected the observation of Wood et al. |2002) in the equatorial plane. Polar starspots may completely change the 
structure of stellar winds in the rapidly rotating, young solar analogs. It is known that the dense slow solar wind is initiated above 
active regions with closed magnetic field lines at low latitudes, while the fast solar wind originates mostly in high-latitude coronal 
holes with opened field lines. The polar appearance of starspots may force the slow stellar wind to flow from the polar regions, so 
almost parallel to the rotation axis. Then this may change the current views about the magnetic braking of stellar rotation in the 
young age of their evolution. 

Stellar magnetic activity has a tremendous influence on the planetary evolution. Coronal mass ejections in young, more active 
stars may lead to the erosion of atmospheres in exoplanets (Khodachenko et al. 120071 Lammer et al. I2007l i. On the other hand, 
the polar spots in young stars could change both the structure of coronal magnetic fields and the formation of prominences, which 
are the actual cause of CMEs. The polar location of prominences may trigger polar CMEs, which propagate significantly out of 
the ecliptic plane. Then the propagation of shocks and magnetic ropes/clouds out of the ecliptic plane may completely change the 
current picture of stellar-planetary interaction. The polar spots and related enhanced plasma flow out of the ecliptic plane during the 
fast-rotation period of the young Sun, or similar stars may lead to a longer stability of atmospheres related to planets that orbit in the 
ecliptic plane, and therefore to more convenient conditions of planetary evolution. Although our findings agree with the mass-loss 
observations of Wood et al. ( 120051 ). clearly more measurements are needed to better define the mass-loss, radiation activity relation 
of cool main-sequence stars, at high activity levels. More observations of G and K-type stars with X-ray surface fluxes greater than 
10^ erg cm - s ' are needed to find out what is happening to stellar plasma outflow during high coronal activity periods. All these 
problems need to be studied in more detail in the future. 

6. Conclusions 

Large-scale magnetic Rossby waves may play a significant role in the dynamics of rapidly rotating stars. Here we have studied the 
shallow-water MHD waves in the tachoclines of rapid rotators. We considered a rotating spherical coordinate system and a toroidal 
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magnetic field in the tachocline. We derived the dispersion relations for high (magnetic Poincare waves) and low (magnetic Rossby 
waves) frequency spectra separately in the approximation of rapid rotation. We found that magnetic Rossby waves tend to locate 
mainly at the polar regions, while magnetic Poincare waves are located at low latitudes. Magnetic Rossby waves may lead to the 
enhancement of magnetic flux at the polar regions and trigger its eruption towards the surface in the form of magnetic flux tubes. 
Then the magnetic tubes may form starspots at the polar region. Magnetic Rossby waves may cause the observed intermediate 
period oscillations in the magnetic activity of rapidly rotating stars, such as young solar proxies. The polar location of prominences 
during the fast-rotating period of the young star may trigger polar CMEs, which propagate significantly out of the ecliptic plane and 
completely change the current picture of the solar/stellar plasma interaction with young planetary environments. 
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